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Abstract 

We study both CP conserving and violating contributions to the decay Kl 
-K^-K^vv. We find that the decay branching ratio is dominated by the CP conserv- 
ing part. In the standard CKM model, we estimate that for mt ~ 174 GeV , the 
branching ratio due to the CP conserving (violating) contributions can be as large 
as 4.4 X 10-13 (^10 X IQ-^^). 



I. INTRODUCTION 



With the prospect of a new generation of ongoing kaon experiments a number 
of rare kaon decays have been suggested to test the Cabibbo-Kobayashi-Maskawa 
(CKM) paradigm: Quarks of different flavor are mixed in the charged weak 
currents by means of an unitary matrix V. However it is sometimes a hard task to 
extract the short-distance contribution, which depends on the CKM matrix, because 
of large theoretical uncertainties in the long-distance contribution to the decays 
1^. To avoid this difficulty, much of recent theoretical as well as experimental 
attention has been on searching for the two modes: —>■ tt^uu and — > tc^uu. 
It is believed that these two decays are free of long-distance and other theoretical 
uncertainties ^ ^f. 

It has been shown that the decay branching ratio of —>■ tt^uu is at the 
level of 10^^'^ 1^, 1^ arising dominated from the short-distance loop contributions 
containing virtual charm and top quarks. This decay is a CP conserving process 
and probably the cleanest one, in the sense of theoretical uncertainties, to study 
the absolute value of the CKM element Vtd- The current experimental limit is 
Br{K+ — > T^^w^expt ^ 5 X 10 ^'^1^ given by the ongoing E787 experiment at BNL. 
It is expected that the experiment will reach the standard model predicted level in 
a few years. On the other hand, the decay Kl ir^uu depending on the imaginary 
part of Vtd is a CP violating process ^ and offers a clear information about the origin 
of CP violation. In the standard model, it is dominated by the Z-penguin and W-box 
loop diagrams with virtual top quark. But there has been no dedicated experimental 
search for this decay yet. Although there are several interesting proposals to study 



this mode at the next round KEK and FNAL experiments |T0[, the experimental 



sensitivities can only be around 10~^, whereas the decay branching ratio in the 
CKM model is at the level of 10^^^ [Q. From an experimental point of view very 
challenging efforts are necessary to perform the experiments. This is because all the 
final state particles are neutral and the only detectable particles are 27's from tt^. 

In this paper, we examine the decay Kl -k^-k^vv. Like the decays of 
TT^uu and Kl -k^vv, we expect that this mode is also a clean one due to the 
absence of photon intermediate states.0 Moreover, in contrast with Kl — > n^uu, it 
contains two charge particles 7r+ and vr" in the final states and could be relatively 
easy to do an experiment Therefore, it should be interesting to give a theoretical 



analysis on this decay to see whether it could be tested experimentally in future kaon 
facilities. 

The paper is organized as follows. In Sect. II, we study the decay rate of 
Kl Ti^n^uu from the short and long distance contributions. We present our 
numerical results in Sect. III. The conclusions are given in Sect. IV. 

II. DECAY RATES 

We start by writing the decay as 

KLiPK) ^ n+{p+)n-ip.),^{k+Mk_) (1) 

where px , p+ , P- , k+ and k_ are the four-momenta of Kl , 7r+ , tt^ , u and z/, 
respectively. Similar to the decays of K~^ n^uu and Kl tt^uu, the short 
distance contributions, arising from the box and penguin loop diagrams with virtual 
charm and top quarks, dominate the decay branching ratio of Kl ti^ti^uu. The 



"We note that the decay of Kl — > tt+tt vi/ is different from that of Kl tt+tt e+e in which 
it is dominated by the long distance due to the photon intermediate states [0 . 
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effective interaction reievant for tlie process is given by ||T 



a 



eff 



y/2 47r sin^ 6 



W i=c,t 



wliere 7]c — 0.71 and r]t — 1 are tlie QCD correction factors |T^, Xi = m^/M^r and 



Cu{Xi) 



3(xi-2), , Xi + 2 

mxi + 



X,; - 1 



Xi - 1 



(3) 



To obtain tlie matrix element, we follow the analysis of Ki^ decays by Pais and 



Treiman |]15|. We define the following combinations of four-momenta: 



P = P++P~, Q = P+-P- 



L = k++k_, N = k: - k_ . 



(4) 



Similar to K14 decays [0, the decay 7r+7r~z/z/ can be kinematically 

parametrized by five variables: = P^, the invariant mass of tt+tt^ pair; Sjy = L^, 
the invariant mass of uu pair; 9^^, the angle between p+ and L as measured in the 
TT+TT^ cm. frame; 6^, the angle between fc+ and P as measured in the cm. frame; 
and (j), the angle between the normals to the tt+tt^ and planes. The ranges of 
the variables are flHl: 



4M2< <Ml, 



< s, <{Mk- ^/^) 
< 9^, < TT, 



< 



< 27r. 



(5) 



respectively. For the hadronic matrix element, we use the standard parametrization: 



-7r-|s7^(l-75)rf|K°> 



K 



FP, + Gg^ + i^e.^p^UPPQ'^ 



(6) 



where the form factors F , G and H can be related by isospin to the corresponding 
form factors in the matrix element of < 7r+7r^|s7^(l — '')^)u\K^ > in decay. 
These form factors have been evaluated in ChPT at order [|l^, [18|. It is found 
that 

Mk 



F = G 



<! TV 

with = 130 MeV. From Eqs. (2) and (6), we obtain the amplitude of the decay 
— n^n^uu for each neutrino flavor as 

V2 47r sm Ow ~^^t 
M^t,,p.UPPQ-] - . (8) 

With the CPT theorem and Kl^K2 + eKi ^ {K'^ - K°)/V2i, we flnd 

A{K,^n^n-uu) = ^--^^ {zG g4-A^AV.(^*: 

V2 4:71 sm Ow Mk ^ 

+A^X\l - p)C,{xi)]] - (9) 

where A = 0.22 is the Cabibbo angle, A , p and t] are the parameters in the Wolfen- 
stein parametrization of the CKM matrix and we have ignored the contribution 
from Ki part because of the smallness of e parameter. In Eq. (9), the terms pro- 
portional to F and H ^ which represent / = s-wave and 1 = 1 p-wave for the tt+tt^ 
system, are CP conserving and that to G, / = 1 p-wave, is CP violating. 

To write the partial decay rate for (1), it is convenient to introduce the following 
combination of kinematic factors and form factors: 

Fi = -iF x\r],C,{x,) + A^\\l- p)C,{xt) 



-a^{P ■ L) COS e^G A'yr]C^{xt) 



F2 = -^^(s^s,)^^ [A2AV.(xt) 
F3 = -a^X{s^s,)^^Uc,{x,) + A^X\l-p)C,{xt) 



(10) 



where 



1/2 



(7, = (l-T^j ' ^ = [(P • i^)' - , P-i^ = -{Ml -S,-S,). (11) 



The differential decay rate is 



X(Tt^I{st^, s^, 9t^, 9^,, (f))dsT^dsi,dcos6T^dcos6i,d(f). (12) 



The dependence of / on 9^ and is given by 

/ = A + ^2 cos 2^,^ + /a sin^ cos 20 + /4 sin 2^,^ cos (!> + h sin cos (f) 

+Iq cos + It sin sin + /g sin 29i, sin + /g sin^ 9^, sin 20 , (13) 

where /i, • • • , /g depend on s^r , Si/, and ^tt- By integrating over the angles, 9i, and 
(f), we obtain 



h - -h 
3 



ii^ir+(ii^2r+iF3r) sin^^. 



(14) 



where we have used the formulas for the form factors /i and I2 given by 



- 4 



i^ir+^(i^2r+iF3nsin^^.' 



(15) 



Combining Eqs. (10)-(14), we get the differential decay rate of — > tt'^^tt for 
three generations of neutrinos as follows 

d^V ( d^V \ ( d^r \ 



ds^^dsi^d cos 9t^ 



ds-j^dsi/d cos ^Ti" J 
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+ 



ds^^dsiid cos 9t^ 



(16) 



CPV 



with 



• (r7ca(xc) + A2A^(1 - p)a(xO)' (17) 



and 

2 



2 



corresponding to the CP conserving and violating contributions, respectively. As 
comparisons, we give the branching ratios for — > n^uu and Kl tt^uu from 
the short distance contributions: 

BriK+^Ti+uu) = -( Br(K+ ^ir'^e+u) 

{vcC,{x,) + A^\\l - p)G,{xt)f + (A^AV- 



2 V27rsin2^VK/ r(i^+) 



•(A^AV.la:*))' , (19) 

where Br{K^ ir^e^u) = 0.048. It is interesting to note that the Kl -k^-k^vv 
decay rate of the CP violating part in Eq. (18) has a similar CKM dependence 
as Kl Ti^uu in Eq. (19) while that of the CP conserving part in Eq. (17) is 
somewhat different from the CP conserving decay of K~^ — * tt+z/i/. 

The long distance contribution can be calculated in the framework of chiral 
perturbation theory. There are three kinds of terms which contribute to the process 
of interest, L'^)^^ 0, reducible anomaly {Lr,a.) and direct anomaly {L^.a) pl| - 

^(2) 



L^f ^ is the weak chiral lagrangian of O(p^) 



^A5=i ^ ^^TrXeD^U^D^U (20) 



where 



is the nonhnear reahzation of the octet meson fields and 



U = expf^^'^A'') (21) 



D^U = d^U -ir^U + iUl^ (22) 



is the covariant derivative with 



= -^zJ-sin^e^Q). (23) 
cos Ow ^ ' 

The overall normalization Gg is determined by the amplitude of K tttt and the 
numerical value is 9 x 10^^ GeV~'^. The matrix Q is the quark charge matrix, 
Q = diag{2/3, —1/3, —1/3), which characterizes the E.M. current coupling of Z. 
The parameter ^ inside the left handed current is the coefficient for the singlet 
current coupling of Z, and it is of unity in the limit of nonet symmetry. Note that we 
have different identification for and than those in The reducible anomaly 
arises from the kind of diagrams starting with a. K — n (or K — if) weak transition 
induced by L^^^, then followed by a tt (or rf) pole and ended by an anomaly vertex 
derived from Lwz.w- The relevant pieces to our calculation in Lwz.w are given by 

'^w.z.w = --^Tre^,o.pL^L''L^f + -^Tre^,^pRm^R^r^ (24) 



where 



= lUd^Ul (25) 



The direct anomaly can be understood as the bosonization of the product of the left 
handed currents arising from L^)^^ and Lw.z.w.- It reads 



^d.a. — 



327r2 

+a2TrXe[U^Fil''U, L^L,] + SasTrX^L^TriF^" + U^F^''U)L, 
+a^TrXeL,Tr{F^'' - U^F^''U)u} (26) 

where F^\ are the field strengths associated with the fields and correspond- 
ingly. The coefficients are in principle of order one and they can be extracted 
from the anomalous radiative decay modes of koan. In terms of the kinetics variables 
defined before, the decay amplitude resulting from long distance effect is given by 

2 

2(3ai — 3a3 — 04) + 2sin^ Oyj(a2 + 204) + £ — k — — 

mj^ - r 

-8i7rV'(l - 2 sin^ e^)(P'' + L'^)}, (27) 



and the corresponding differential decay rate is then given by 

+al sin^ OwS-j^Sj, [2(3ai — 3a3 — 04) 
+2sm^e^{a2 + 2a^)+im\/{m]^-ml)^Y (28) 
III. NUMERICAL RESULTS 

The validity of relating rrit to the decay rate depends upon the negfigibility of 
long distance contribution. Therefore it is important to learn the branching ratio 
arising from the long distance effect. Due to the absence of mt in the amplitude, 



9 



the decay rate of long distance contribution is relatively suppressed by at least two 
orders. Numerically we find 



-20 



Br{Lfi=^) = 5.0 X 10-^^ (29) 



Br{Ld.a. + Lr.a) ~ 4.7 X 10 

- 

-(2) 



As we shall see below it is safe to ignore the long distance effect and we shall 
concentrate on the short distance effect only in the following analysis of decay rate. 

To estimate the CP conserving and violating decay rates in (17) and (18), we need 
to find out the allowed values for the CKM parameters A , p and rj, constrained by 
the experimental measurements such as e, the CP violation parameter in ^ vrvr; 
Xd, the — B^ mixing; and the ratios | Vcf,/\4^s| and |\4b/Vcfe| of the CKM elements. 



We use the same fitting procedure and the necessary equations in Refs. fg, 

the fits, we take the updated values iKbl = 0.041±0.005 and IK^/Ktl = 0.080±0.025 

and /b = 200 ± 50 MeV. 

Integrating over all the variables in Eqs. (17) and (18), we can examine the 
decay rates for both CP conserving and violating parts which depend on the top 
quark mass and the CKM parameters. In Figures la and lb, we plot the branching 
ratios of CP conserving and violating contributions to tt^tt^uu as a function 

of the top quark mass, showing the lower and higher values allowed at 90% C.L., 
where Br{KL -k^-k^vv) = T{Kl 'k^'k^vi')/T{Kl all). We also show the 
corresponding decay branching ratios of n^uu and —>■ tt^uu in Figures 

Ic and Id by using Eq. (19), respectively. From the figure, we see that the CP 
conserving part of the branching ratio is much larger than that of CP violating one. 
Clearly, measuring the decay rate will not give us information on the CP violation. 
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For 150 <mt< 200 GeV, we find 

1.1 X 10-^3 < Br{KL n+n-up)cpc < 5.0 x 10"^^ , 

5.0 X 10-^^ < Br{KL n+n-uP)cpv < 1-1 x 10"^^ . (30) 

We now study the differential decay spectrum in terms of {9^^) by integrating 
over Si, and 9^, (s^r and s^) in Eqs. (17) and (18) to see whether we would distinguish 
the CP conserving and violating parts. We define the normalized invariant mass of 
7r+7r~ as X — s^/Ml-. To illustrate the shapes of the spectra between the CP 
conserving and violating cases, we choose rrit ~ 160 GeV and the CKM parameters 
A ~ 1 , p ~ —0.2 and rj ~ 0.4. We plot the differential branching ratios dBr{KL — > 
TT~^TT~i'u)/dx^ vs and dBr{KL —>■ tt^t^'vv) / dcos9-rf vs cos6'7r in Figures 2 and 
3, respectively. As shown in Figure 2, the CP conserving and violating spectra of 
dBr^Ki — > 7r''"7r~z/z/)/(ia;^ have similar shapes and are dominated by small values of 
Stt- However, in Figure 3, as expected, [dBr{KL 7r~^ tt' i/i?) /d cos 9Tr](~,py becomes 
maximum when 9^^ is close to or tt and to minimal when it reaches 7r/2 whereas 
[dBr{KL t:^ TT' vv) /d cos 9 ■j^]qpq does the opposite way. Unfortunately, the values 
of the CP violating one around = 0, tt may be still too small to be tested. 

IV. CONCLUSIONS 

We have studied both short and long distance contributions to the decay of 
Kl TT~^TT~i'u. We have demonstrated that the long distance effect to the decay rate 
is neghgible small. We have shown that the branching ratio of the decay is dominated 
by the CP conserving part. With the updated CKM parameters, we find that the 
decay branching ratio is predicted to be (1—5) x 10~^^ for rrit < 200 GeV, which could 
be accessible to experiments at future kaon facihties. The CP violating contribution 

11 



to the branching ratio seems impossible to be measured in experiments. However, 
it is, in principle, to distinguish the CP conserving and violating contributions by 
measuring the spectra of the angular dependence of the differential decay rates. 
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Figure Captions 

Figure 1: Allowed branching ratios for (a) CP conserving contribution to — > 

7r"'"7r~i/P; (b) CP violating contribution to tt'^tt'i/P; (c) — > 
TT+i/P; and (d) Kl n^uP as functions of nit at 90% C.L. 

Figure 2: The differential decay spectrum of dr{KL — > n'^7r~i'P)/dx^ as a func- 
tion of x'2 = y^/MK with mt = 160 GeV , ^ 1.0 , p 0.2 and 

77 ~ 0.4. 

Figure 3: The differential decay spectrum of dr{KL 7r'^7r~uP)/dcos9-!^ as a 
function of cos6'7r. Legend is the same as in Figure 2. 
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